Evolution of the persistent spin helix in the presence of Hartree-Fock fields 
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We derive a spin diffusion equation for a spin-orbit coupled two-dimensional electron gas including 
the Hartree-Fock field resulting from 1st order electron-electron interactions. We find that the life- 
time of the persistent spin helix, which emerges for equal linear Rashba- and Dresselhaus spin-orbit 
interactions, can be enhanced considerably for large initial spin polarizations due to the Hartree- 
Fock field. The reason is a reduction of the symmetry-breaking cubic Dresselhaus scattering rate by 
the Hartree-Fock field. Also higher harmonics are generated and the polarization of the persistent 
spin helix rotates out of the (S y ,S z )-pl&ne acquiring a finite ^-component. This effect becomes 
more pronounced, when the cubic Dresselhaus spin-orbit interaction is large. 



I. INTRODUCTION 

In recent years the spin dynamics in semiconductors 
has attracted great interest, because integrating the spin 
degree of freedom into semiconducting devices opens up 
new possibilities of realizing spintronics or quantum com- 
putational devices.^ For these potential applications long 
spin lifetimes are highly desirable. In this respect the per- 
sistent spin helix (PSH), 2 3 a spin-density wave state with 
ideally infinite lifetime, is of great interest. This state 
exists in two-dimensional electron systems with Rashba 
and linear Dresselhaus spin orbit interactions of equal 
magnitude due to a SU(2) symmetry of the underlying 
HamiltonianP On a less abstract level, the PSH state 
emerges because for Rashba and linear Dresselhaus spin 
orbit interactions of equal magnitude the resulting spin 
orbit interaction points into the same direction for all 
electronic momentum states. The magnitude of the spin- 
orbit interaction grows linearly with the projection of the 
momentum on the PSH wave-vector, such, that all spins 
rotate by the same angles when they diffuse along the 
PSH direction. 

In real systems the lifetime of the PSH is no longer in- 
finite due to the presence of additional terms which break 
the SU(2) symmetry of the Hamiltonian. It has been sug- 
gested, that the presence of cubic Dresselhaus spin orbit 
interaction is the dominant symmetry-breaking mecha- 
nism leading to a finite lifetime of the PSH statel^liSl 
Experimentally, the PSH state has been observed in 
GaAs/AlGaAs quantum wellsj^' where the cubic Dres- 
selhaus spin orbit interaction is quite small, as well as in 
(In,Ga)As quantum wellspl where the cubic Dresselhaus 
spin orbit interaction is much larger. 

The observed temperature dependence of the life- 
time of the PSH state^ suggests, that also electron- 
electron i nterac tions strongly affect the lifetime of the 
PSH stateP^LUJln our previous worle^we found that the 
consideration of electron-electron interaction in combina- 
tion with cubic Dresselhaus spin orbit interaction as the 
symmetry breaking mechanism can fairly well account 
for the observed temperature dependence of the PSH life- 
time. In this previous treatment we have only considered 



electron-electron collisions, which arise in 2nd order in 
the electron-electron interaction but we have neglected 
the Hartree-Fock term, which arises in 1st order. This 
approximation is valid, when the initial spin polarization 
is small and seems to be appropriate for the present ex- 
perimental observations of the PSH stateP On the other 
hand, it is experimentally also possible to realize large 
initial spin polarizations, 1 ^ where the Hartree-Fock term, 
which acts like an effective magnetic field pointing along 
the local spin polarization, is known to enhance the spin 
lifetime considerablyP^HH Therefore it is the purpose of 
the present paper to analyze the effect of the Hartree- 
Fock term on the PSH in the diffusive D'yakonov Perel' 
regime. In particular we want to answer the following 
questions: (i) what is the effect the Hartree-Fock term 
on the lifetime of the PSH-state and (ii) does it modify 
the pattern of the spin-polarization of the PSH-state? 

This article is organized as follows: In Sec. |III| we 
present the spin diffusion equation for a two-dimensional 
electron gas with Rashba- and Dresselhaus spin-orbit in- 
teractions including the Hartree-Fock field, which arises 
from 1st order electron-electron interaction. In Sec. IIVI 
we show how the presence of the Hartree-Fock term en- 
hances the lifetime of the PSH state. Furthermore we 
demonstrate that the presence of the Hartree-Fock term 
also generates higher harmonics in the spin polarization. 
In particular, we show that the 2nd harmonic corre- 
sponds to an S^-component and thus rotates the spin- 
polarization out of the (S y , S z )- plane. The appendices 
contain details about the derivation of the Hartree-Fock 
term (appendix [A]) and of the spin diffusion equation 
(appendix [Bj. 



II. MODEL 

We describe conduction band electrons in a (001) 
grown quantum well with Rashba- and Dresselhaus spin- 
orbit interactions, potential impurities and electron- 
electron interactions. For simplicity we focus on the 
limit of zero temperature, where the 2nd order electron- 
electron scattering vanishes but the 1st order contribu- 
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tion due to electron-electron interaction, i.e., the Hartree- 
Fock term, is finite. Although our treatment is strictly 
valid only at zero temperature, we expect, that qualita- 
tively similar results hold also for small but finite tem- 
peratures, because the temperature dependence of the 
Hartree-Fock term is weak and the presence of 2nd or- 
der electron-electron interactions does not qualitatively 
affect the PSH stateP 

We consider the following Hamiltonian (here, and in 
the following use H = 1): 



H — H + H imp + H e _ e , 



(1) 



where the first term represents a two-dimensional elec- 
tron gas (2DEG) with a quadratic dispersion and intrin- 
sic spin-orbit interaction 



H 



^ ks / Hqs's tpka 

s,s';k 



and the 2x2 matrix in spin space 

H = e k + b(k) • <r. 



(2) 



(3) 



The ip^ s (ipks) are creation (annihilation) operators for 
electrons with momentum k and spin projection s. 
Within the standard envelope function approximation^ 

k 2 

one finds ez, = where m is the effective mass. The 
vector of Pauli matrices is denoted by <j and the in-plane 
spin-orbit field 

b(k) = b fl (k)+b D (k)+b cD (k) (4) 

contains Rashba spin-orbit interaction 



bij(k) = av F ( -k x 




(5) 



Finally, the Hamiltonian ([I]) contains electron-electron 
interactions, 
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H °- C= 2V ^ yk3 - k ^ kl ' k2 ^k 4S2 V'l3 Sl V'k lSl i/'k 2S2 (8) 

ki...k 4 

S1.S2 

with a Thomas-Fermi screened Coulomb potential of the 
form Vk 3 ,k 4 ,ki,k 2 ~ «(l k 3 - ki|)(5 kl+k2 _ k3 _ k4i0 where 



v(q) 



— 27r . and e 



1 



q a* 



with 



4 7T £q S r 

m e 2 



being the effective Bohr radius. 

Following the procedure of our previous treatment of 
the PSH-statePwe use the Nonequilibrium statistical op- 
erator methocP^to derive kinetic equations for the charge 
and spin components of the Wigner-transformed density 
matrix 



/5 k (x, t) = n k (x, t) + s k (x, t) ■ a. 



(9) 



Since, to zeroth order in b/Ep, the kinetic equations for 
spin and charge decouple we only need to consider the 
spin component and find for the spin density s k : 

d t s k + v • 9 x s k + 2 s k x b(k) = J™ p + J k c(1) (10) 

with Vi = ki/m, where the index i = x,y labels the in- 
plane spatial directions. For simplicity, we focus here 
on the limit of zero temperature, where the 2nd or- 
der electron-electron collision integral vanishes and we 
keep only the first order electron-electron collision inte- 
gralJ k ° (1) . As we show in appendix [a] J k °^ 1 ' 1 can be cast 
into an effective magnetic field, i.e., a Hartree-Fock term, 
arising from the mean polarization of the PSH-state itself 



0(1) _ 



-2s k x 



dq 

^n<?)s k +q- 



(27T) 



(11) 



The impurity collision integral in Eq. (JlOj) is given by: 
J k" P = - E ^kk'<5(e k - e k ')(sk - a*) , (I 2 ) 



as well as linear and cubic Dresselhaus spin-orbit inter- 
actions 



ky 

bc(k) = (3v F | k x 




b cD (k.) = 7«f- 




(6) 



Here, vf is the Fermi velocity, the angle 9 gives the di- 
rection of k with respect to the x axis, which we choose 
to point along the (llO)-crystal direction. 

Furthermore, we have included in Eq. ([TJ) scalar 
electron-impurity interactions, 



Simp = y E Wc's^k'kV'ks, 



(7) 



s;k,k' 



arising from scattering of electrons at potential impuri- 
ties V lmp (x) = X]i' y ( x — R*)) where v(x) denotes the 
potential of each single impurity and {R, } are the impu- 
rity positions, eventually to be averaged over. 



where the transition rate is given by 

2irrii\v (k' — k) | 2 with the impurity concentration rij. 



III. SPIN-DIFFUSION EQUATION INCLUDING 
HARTREE-FOCK FIELDS 



In the diffusive strong scattering D'yakonov-Perel' 
regime b{kp)r <C 1 (where r is the momentum relax- 
ation time) one can derive a diffusion equation for the 
spin density of the form (see appendix [B] for details) 



d t S = DS + H : + H 3 



(13) 



We focus for simplicity on equal Rashba and Dresselhaus 
spin-orbit interactions a = /3, i.e., the usual condition for 
the PSH state. In this case we find for the matrix D 



D 



fDdl 



-7cd 

D{dl-ql)~% d 2q Dd y 
-2q Dd y b{d 2 y -ql)-2^ cd/ 



(14) 
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where go = ^rnvpa is the PSH wave- vector. This corre- 
sponds to the usual spin diffusion matrix in the absence 
of Hartree-Fock fields but with a renormalizcd diffusion 
constant D and a renormalized cubic Dresselhaus scat- 
tering rate 7 c d, given by 



D = 



D 



1 + (2S 1 r) 2 ' 



7cd = 



7cd 



1 + (2B 3 r)2 



(15) 



where D = \v f t is the ordinary diffusion constant and 
7cd = ^Vp-f 2 kpT is the cubic Dresselhaus scattering rate. 
We also introduced the effective Hartrcc-Fock fields 



Bi,3 — Xi.sS, 



(16) 



which are the Hartree-Fock fields for the winding number 
one (s±i) and winding number three (s±s) components 
of the spin density 



s k = J2 Mk)e lW . 



(17) 



;=o.±i,±3 



For xi,3 we find ( see appendix |A| and appendix [B]) 

(18) 

where 6> q k is the angle between k and q. Assuming a 
Fermi temperature of Tp — AOOK, an effective mass of 
"m c s = 0.067m e and a dielectric constant of e r = 12.9 
(for GaAs) one finds xi = Xk F ,\ = 34.7cm 2 /s and x 3 = 
Xfc F , 3 =43.2cm 2 /s at T = 0. 



The two other terms in (131 read in the limit a = 



Hi = 2D{B lT x (figS) + ql(B lT ■ e x )(S 



x e x 



(19) 



+ q d y [BiT x (S x e x )] + q e x x (d y S x Bit) 
+ 2q Q (B 1 T-d y S)(B 1 T x e x ) 
+ 2d v [B 1 T(B 1 T-d y S)]}, 
H 3 = - 2% d {B Zx r, B 3y r, 0) x S. (20) 

Due to the presence of the terms H13 the spin diffu- 



sion equation ( 13 ) becomes a nonlinear partial differential 
equation, because the Hartree-Fock fields 3 depend on 
the solution for the spin density S (see Eq. (fl6|). The 



exact solution can thus only be found numerically. Note, 
that for the derivation of the spin diffusion equation only 
b(kp)r <C 1 was required but no restriction applies with 
respect to the magnitude of the Hartree-Fock fields -Bi j3 , 
since those fields have winding number zero. 

In our previous treatment)^ of the PSH state in the 
presence of cubic Dresselhaus spin-orbit interaction we 
have found that the PSH state survives in the presence 
of cubic Dresselhaus spin-orbit interaction but the am- 
plitudes of the spin polarization of the PSH-state acquire 
a small ellipticity and a finite lifetime. The PSH eigen- 



mode of ( 14 1 in the absence of the Hartree-Fock field (i.e. 



with Bi y3 —> 0) reads 
S x = 

S y = -So sin(<7 y)e- t/TpSH 



(21) 



*: = ■% I \i 1 + ^ _ cos(goy)e- t/rpSH 



So 1- 



7cd 



cos(g y)e 



"*/ T'PSH 



where the lifetime tpsh is given by 

^ = 5r + § 7od -i^S«§ 7od . (22) 

Here, we have denoted the D'yakonov-Perel' relaxation 
rate by T — iq^D and the approximate expressions 
for S z and tpsh are valid to lowest order in 7cd/r oc 
[6 c d(^f) /bR(kp)} 2 , which is basically the square of the 
ratio of the cubic Dresselhaus spin orbit interaction over 
the linear spin orbit interactions. 

The purpose of the present paper is to investigate how 
the pattern and the lifetime tpsh of the PSH state are 
modified in the presence of the Hartree-Fock field aris- 
ing from first-order electron-electron interactions. In 
Sec. |IV A"| we present our results concerning the lifetime, 
based on numerical solution of the diffusion equation 
Then, in Sec. |IVB we present a more analyt- 
We will simplify the prob- 



(131 



(13) 



ical analysis of Eq. 
lem and investigate the diffusion equation assuming (i) 
the limit of small cubic Dresselhaus spin-orbit interac- 
tion, 7cd/r <C 1, and (ii) a time-independent helical 
magnetic field instead of the self-consistent Hartree-Fock 
field. These assumptions allow for an analytical ap- 
proach, making use of the smallness of the parameter 
7cd/r. Although corresponding to a qualitatively differ- 
ent situation, this approach provides some insight in the 
general structure of the solutions of Eq. (13 \ . Finally, in 



Sec. IV C we present full numerical solutions of Eq. ( 13 ) 



now focusing on the spatial pattern of the PSH for dif- 
ferent initial spin polarizations and for small and large 
cubic Dresselhaus spin-orbit interaction. 



IV. MODIFICATIONS OF THE PSH STATE 



Lifetime enhancement of the PSH state 



Our main finding is, similar to previous workp^UHl 
that the PSH lifetime is enhanced considerably due to 
the Hartree-Fock field arising from first-order electron- 
electron interactions. The intuitive explanation is simply 
that the Hartree-Fock field is parallel to the spin polar- 
ization of the helix state and thus strengthens the spin 
polarization. Since, experimentally this can e.g. be real- 
ized by increasing the initial spin polarization, i.e., the 
amplitude of the PSH state, we show in Fig. [T] the time 
evolution of the S z spin polarization for initial spin polar- 
izations of the PSH state with 10% (blue) and 20% (red) 



4 



B. Analytical approach for 7cd/T <C 1 



without Hartree Fock 
10% spin polarization 
20% spin polarization 




PSH 



FIG. 1. Time dependence of the ^-amplitude of the PSH- 
state for different initial spin polarizations (solid lines). For 
comparison the time evolution without Hartree-Fock field is 
shown in black. Dashed lines indicate, how the time depen- 
dence of the Sz-amplitude of the PSH-state would look like, if 
the Hartree-Fock field remained constant at its initial value. 



using the parameters of the quantum well of Ref. |U a 
Fermi temperature of Tp — AQOK, an effective electronic 
mass of ?7i ff = 0.067m e , linear Rashba and Dresselhaus 
spin orbit couplings of a = j3 = 0.0013 and a cubic Dres- 
selhaus spin-orbit interaction of jvp = 5.0eVA 3 as well 
as a momentum relaxation time of r = lps. In the sim- 
ulations we always assume that at time t = the elec- 
tronic spin polarization corresponds to the PSH state of 
Eq. (21 ). We then solve for the time-dependent polariza- 



tion S[i) with the Hartree-Fock held switched on at time 
t = 0. For comparison the time evolution of the S z spin 
polarization of the PSH state without Hartree-Fock field 
is also shown (in black) . Obviously, for realistic quantum 
well parameters the lifetime of the PSH state can easily 
be enhanced by factors of 10 to 20 when the initial spin 
polarization is in the percentage range. 

The main cause for the lifetime enhancement of the 
PSH-state can be traced back to the reduction of the ef- 



fective cubic Dresselhaus scattering rate 7 c( j in Eq. ( 15 ) 
in the presence of a Hartree-Fock field. Since it is the 
SU(2)-symmetry breaking cubic Dresselhaus spin-orbit 
interaction, that renders the lifetime of the PSH state 
finite in the first place, a reduction of the effective cu- 
bic Dresselhaus scattering rate considerably enhances the 
lifetime of the PSH state. Because the Hartree-Fock field 
is proportional to the spin polarization itself, its mag- 
nitude and therefore its ability to suppress the effective 
cubic Dresselhaus scattering rate 7 C( j decreases with de- 
creasing spin polarization. As a consequence the spin 
polarization (solid lines) in Fig. [I] decays faster than the 
exponential decay (dashed lines) expected for a constant 
field. 



Whereas the full nonlinear solution for the spin polar- 
ization in the presence of a self-consistent Hartree-Fock 
field can only be obtained numerically, we can simplify 
the problem and work with a static magnetic field in- 
stead. We take as an ansatz for the PSH state 





>i>sii = S ( -singoy 

(1 - x) cos q y. 



(23) 



where x parametrizes the ellipticity of the spin polariza- 
tion due to the presence of cubic Dresselhaus spin-orbit 



interaction (see Eq. (21)). We assume the fields B13 to 



be static and proportional to this polarization profile, 

Spsh 



Bi 3 = B\ 



'PSH I 



(24) 



diffusion equation (13 1, we find that with 

1 



If we plug the ansatz (|23|) and the helds |24|) into the 

(25) 



7cd 



1 



AqiDl + ^Bir) 2 



Ted 

f l + i(Sir) 2 ' 



the instanteneous time-evolution of the PSH state reads 
to lowest order in 7 c d/r 



d t S 



't^PSH 



»PSH 

•?PSH 



&2x 

Tlx 



&3yz 
T 3yz 




3f PSH l + i(5iT) 2 





- sin(3g 2/) 
cos(3<7o2/) 



(26) 

(27) 
(28) 

(29) 



This equation thus describes two effects: (i) The state 
Spsh taken as ansatz tends to relax towards S = with 
an instantaneous relaxation time fpsHj which is enhanced 
in comparison to the bare tpsh due to the reduction of 
the cubic Dresselhaus scattering rate in the presence of a 
helical magnetic field, (ii) The polarization starts to de- 
velop higher harmonics (with wave vectors 2qo and Sqo), 
including an x-component which was not present before 
(cf. Eq. m). 



As we will show below, the magnitudes of the higher 
harmonics are smaller by factors 7cd/r (see Eq. (33)). 
This means that in the limit of 7cd/r < 1 wc can ne- 
glect the contributions of S2 X and S^ yz in (261 to first 



approximation, and Eq. ( 26 ) becomes valid at all times 
(at least for constant magnetic fields -B1.3). In this case 



5 



one can solve the equation easily to find the simple ex- 
ponential decay 



Spsh(^) — SpsH e 



-t/fpSH 



(30) 



This exponential solution for constant fields is indicated 
by dashed lines in Fig. [T] For the actual time-dependent 
Hartree-Fock field (solid lines in Fig. [T]) the lifetime is of 
course shorter than tpsh of Eq. ( p7| : as soon as the spin 
polarization of the PSH state starts to decay, the fields 
B1.3 decrease proportionally, which results in an increase 
of the decay rate 1/fpsH of the PSH state. 

Let us now turn to the generation of higher harmonics 
in ( 26 ) , which is due to the nonlinearity of the spin diffu- 



sion equation. In particular a second harmonic sin(2g y)- 
term appears in the S x component (see Eq. (28)) and a 
third harmonic sin(3(7o2/)-, cos(3go z )-term appears in the 



S y - and ^-components (see Eq. (29 1). Of course, this is 



not the end of the story, since, in order to obtain the full 
solution, the second and third harmonics would have to 
be included into the initial expression leading to a suc- 
cessive generation of all higher harmonics. 

The most surprising and at first counterintuitive find- 
ing is the generation of a finite S^-component of the spin 
helix, i.e., the observation that a Hartree-Fock field par- 
allel to the local spin polarization rotates the spin po- 
larization of the PSH out of the yz-plane by creating 
an additional out-of-plane ^-component. In order to 
estimate the magnitude of this new second harmonic S x - 
component of the PSH-state, we modify the ansatz of 



the initial PSH state of Eq. (23 1 to include the second 
harmonic ^-component as 



s x sin 2q y 
Spsmx = S | - sin q y 

(1 - x) cos q y. 



(31) 



Plugging this ansatz into the spin diffusion equation ( 13 ) 



one finds for the ir-component of SpsH2x to leading order 
in 7cd/r 



d t (S s x ) = So-y 1 ( 



BoT- 



1 



Bit 



(32) 



We see that the second harmonic has a fast relaxation 
rate T compared to the time scale l/j c d of the other 
terms. Let us solve (32) for s x (0) — 0, i.e., starting 
from the PSH ansatz of Eq. (23). Assuming that the 



spin polarization of the PSH state So remains constant 
in time we find 



Ted 
2f 



B 3 t 



B lT 



1 



\ (Sir) 2 



4ft 



(33) 



From this simple consideration we see that the fast re- 
laxation rate of s x has two effects: (i) It makes the equi- 
librium amplitude of the S^- component very small for 
small cubic Dresselhaus spin-orbit interaction, since it is 
proportional to 7 c d/r. (ii) It causes the ^-component 




FIG. 2. For equal amplitudes S(y + dy) = S(y — dy) 
the contribution of the winding number one spins to the S x - 
component diffusing from y + dy and from y — dy to the point 
y would cancel, because they rotate around the same local 
Hartree Fock field B(y) but correspond to ±fc-states and thus 
are oriented oppositely. Due to the elliptical profile of the 
PSH amplitude the amplitudes S(y + dy) and S(y — dy) differ 
and thus a finite ^-contribution emerges. 



to rise from zero to its maximum value within a short 
switch on period given by 1/r <C tpsh- 

Can we understand the physical origin of this sec- 
ond harmonic S^-componcnt? Normally, one would ex- 
pect that a parallel magnetic field simply strengthens the 
parallel spin orientation. In the presence of spin-orbit 
fields, however, the spin density consists of parts with 
different winding numbers, where only winding number 
zero contributes to the local spin polarization while the 
other winding numbers (one and three) average out to 
zero. Since the anisotropic components of the spin den- 
sity (winding numbers one and three) are not necessarily 
parallel to the local spin density and thus not necessarily 
parallel to the Hartree-Fock field, they can modify the 
pattern of the PSH state by rotations around the local 
Hartree-Fock field. 

The expression for the second harmonic S^-component 
in the presence of a constant Hartree-Fock field (see 
Eq. (33)) suggests that both winding number one and 



winding number three parts of the spin density contribute 
to the formation of the S^-component. The winding 
number one part contributes only via the diffusion term. 
At one fixed position y electronic spins arrive via diffusion 
from y+Sy with amplitude S(y+Sy) and from y—5y with 
amplitude S(y — Sy). While diffusing to position y both 
spins arriving from y±Sy rotate around the local Hartree- 
Fock field, which is parallel to S(y) and thus acquire a 
finite ^-component, since S is polarized in the yz-plane 
(See Fig.[2]). Since the angle between S(y+Sy) and S(y) is 
equal to the angle between S(y) and S(y — Sy) the result- 
ing ^-components are of opposite sign an would cancel 
for equal magnitudes S(y + Sy) = S(y — Sy), i.e., for a 
PSH state with equal S y and S z amplitudes. However, in 
the presence of cubic Dresselhaus spin orbit interaction 
the initial PSH profile is elliptical with the S y and S z 
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FIG. 3. Time dependence of the S^-amplitude, which 
arises solely due to the Hartree Fock term, at first maximum 
qoy = 7r/2 (left) and spatial and temporal evolution of the S x - 
amplitude of the PSH state (right) for an initial spin polar- 
ization of 20% and with linear spin-orbit interactions reduced 
by a factor of 0.2. (7cd/T = 0.067 and r PSH = 412ps). 



FIG. 4. Time dependence of the ^-amplitude at first max- 
imum qoy = tv/2 (left) and spatial and temporal evolution of 
the ^-amplitude (right) of the PSH state for an initial spin 
polarization of 20% and with linear spin-orbit interactions re- 
duced by a factor of 0.02. (7 c d/r = 6.8 and tpsh = 324ps). 



amplitudes differing by approximately 7cd/r. Thus the 
amplitudes of the electronic spins arriving from y ± 5y 
differ due to the elliptical PSH profile leading to the for- 
mation of a finite ^-component (see Fig. [2]). Only when 
the local spin polarization is parallel to one of symmetry 
axes of the ellipse, i.e., for S y — or S z — the ampli- 
tude of the electronic spins arriving from y±5y are equal 
and thus S x — at these points. Since this situation oc- 
curs for qoy — 0, 7r/2, it, 3tt/2 this explains the oscillation 
of S x with double wave vector 2qo. 

The winding number three part of the spin distri- 
bution function contributes to the formation of an S x - 
component via combined rotations around the local 
Hartree-Fock field and the component of the cubic Dres- 
selhaus spin orbit field parallel to e y . Again a finite 
contribution to S x arises only when the local spin po- 
larization is neither parallel nor perpendicular to e y (the 
direction of the relevant component of the cubic Dres- 
selhaus spin orbit field) thus leading to an oscillation of 
the winding number three contribution to S x with dou- 
ble wave vector 2qg similar to the winding number one 
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case 1 



C. Modification of PSH pattern 

Since the magnitude of the S^-component as well as 
the other higher harmonics are of higher order in 7 c d/r, 
they will become significant for systems with large cu- 
bic Dresselhaus spin orbit interaction. Therefore it is 
interesting to investigate how the spatial and temporal 
evolution of the PSH state is modified when cubic Dres- 
selhaus spin orbit interaction is large 7 c d « T. Since the 
magnitude of the cubic Dresselhaus spin orbit interaction 
is given by the crystal symmetry of GaAs, whereas the 
magnitudes of the linear Rashba and Dresselhaus spin 
orbit interactions a and j3 can be varied by changing the 
doping asymmetry and the width of the quantum well, we 
could access the 7 C( j ~ T-regime by reducing the magni- 
tude of the linear spin orbit interactions while keeping the 
magnitude of the cubic Dresselhaus spin-orbit interaction 
fixed. Since the lifetime of the PSH depends mainly on 




FIG. 5. Spatial and temporal evolution of the .S^-component 
(left) and the S'y-component (right) of the PSH state for the 
same parameters as used in Fig. [4] 



the magnitude of cubic Dresselhaus spin-orbit interaction 



(see Eq. (21), where one finds r P g H w §7 c d for 7 C( j <C T 
and Tpg H w 7 C( j for 7 c d 3> T), this does not significantly 
affect the absolute lifetime of the PSH state. Alterna- 
tively, the regime of large cubic Dresselhaus spin-orbit 
interaction could also be reached by using an (In,Ga)As 
quantum well instead, or by increasing the Fermi energy. 

In Fig. [3] we show the temporal and spatial evolution of 
the ^-component of the PSH state, which arises solely 
due to the Hartree-Fock term, for only slightly enhanced 
cubic Dresselhaus spin orbit interaction. This situation 
has been reached by a reduction of the linear spin orbit 
interactions by a factor of 0.2, which results in a ratio of 
7cd/P = 0.067. As expected from Eq. (33) the magnitude 
of the S^-component is quite small and its spatial vari- 
ation is characterized by a sin(2goy)-dependence. The 
initial switch on process now exhibits an oscillatory be- 
havior and, contrary to the S y - and ^-components, the 
magnitude of the S^-component saturates before it drops 
to zero quite abruptly. 

As an example for large cubic Dresselhaus spin orbit 
interaction we show in Fig. [4] the spatial and temporal 
evolution of the PSH state for linear spin-orbit interac- 
tions reduced by a factor of 0.02, which results in a ratio 
°f 7cd/r = 6.8. The magnitude of the ^-component is 
now quite sizable, reaching almost 8% of the PSH ampli- 
tude at intermediate times. Obviously, the spatial vari- 
ation is no longer characterized by a simple sm(2qoy)- 
dependence but even higher harmonics start to play a 
prominent role. For comparison the spatial and tempo- 
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ral evolution of the S y - and ^-components are shown in 
Fig. 5] Due to the large ratio of 7cd/r = 6.8 also the 
amplitude of the S z -component is strongly reduced and 
the presence of higher harmonics is clearly visisble. 



V. CONCLUSIONS 

Based on a spin-coherent Boltzmann type approach we 
have derived a semiclassical spin diffusion equation for a 
two-dimensional electron gas with Rashba- and Dressel- 
haus spin-orbit interactions including the Hartree-Fock 
field, which arises in first order in the electron-electron 
interactions. In agreement with previous treatments"^ 
we find that the Hartree-Fock term considerably en- 
hances the lifetime of the persistent spin helix state, when 
the initial spin polarization is in the percentage range. 
Within our approach we are able to assign the origin of 
the lifetime enhancement to the suppression of the ef- 
fective cubic Dresselhaus scattering rate by the Hartree 
Fock field. 

Surprisingly, not only the lifetime but also the pattern 
of the spin polarization is modified in the presence of a 
Hartree Fock field. Due to the nonlinearity of the prob- 
lem higher harmonics in the spin polarization are gener- 
ated. The second harmonic corresponds to a spin polar- 
ization in S^-direction and thus rotates the spin polariza- 
tion of the persistent spin helix state out of the (S y , S z )~ 
plane. These modifications of the spin pattern of the 
persistent spin helix become significant when the cubic 
Dresselhaus spin-orbit interaction is of similar magnitude 
as the linear spin-orbit interactions. Since the magnitude 
of the cubic Dresselhaus spin-orbit interaction is fixed by 
crystal symmetry, we have suggested that this regime 
could be experimentally reached by reducing the magni- 
tude of the linear spin-orbit interactions via appropriate 
design of the quantum wells. 
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2nd order electron-electron interaction vanishes we con- 
sider the kinetic equation only up to first order in V, 
yielding: 



with the first order mean field term: 

4){t) = i{[v,4i, l \) t TeX . 



(Al) 



(A2) 



We evaluate the mean field term <[A2j) for the electron- 
electron interactions of Eq. Q. Evaluation of the com- 
mutator in Eq. ( A2 1 yields in this notation 



JlfU{[H c . c ^\4i]) (A3) 



2'3'23 



Upon Wick decomposition of the averages 

(vM^) = -(4 v 8 ) (4 + (^4) (4h 

= -^k 1 k3 < 5k 2 k4/siS3(kl)/s 2 s 4 (k2) 

+<5k 1 k 4 <5k 2 k 3 f^ s i ( k i)/s 2 s 3 ( k 2) (A4) 

and by exploiting the Kronecker symbols we obtain a 
mean field term that is diagonal in momentum but re- 
mains a matrix in spin space 

C (1) ( k )=^E«(i k '- k i)[/( k ')/( k ) - /( k )/( k ')W ■ 



(A5) 



By decomposing the spin space matrices according to 



/(k) = n k + er-s k . 



3(1) _ 



J 



3(1) 



j; 



e(l) 



(A6) 



one finds that the mean field term entering on the right- 
hand side of the kinetic equation for the spin density 
(Eq. |l0|) becomes: 



3(1) 



-2sk x 



v(q)s k+ci . 



(2tt 



(A7) 
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Appendix A: Hartree Fock term 

For a general nonequilibrium problem with the Hamil- 
tonian H = Hq + V containing an exactly solvable 
single-particle part Hq and an interaction V, Zubarev's 
Nonequilibrium statistical operator formalism 16 allows 
to derive a closed set of equations that describe the 
irreversible temporal evolution of the density matrix 
(tpj/tpi} 1 = T r [/°(£) V'i'V'z]- Here I, I' are possibly composed 
indices, e.g., for momentum and spin. For details of the 
derivation, see Ref. Since for zero temperature the 



Appendix B: Spin diffusion equation 

We follow the approach used in our previous work^ 
to set up a diffusion equation for the spin density Sk 
valid in the diffusive D'yakonov-Perel' regime, i.e., in the 
regime 6i?r -C 1 where due to strong scattering the spin 
polarization is stabilized. In order to solve equation ( 10 ) 



for the spin density we expand the spin density into k- 
space winding numbers 



s k = ^2 s i{k)e' 

Z=0,±l,±3 



10 



(Bl) 



where we keep besides the isotropic component Sq = 
— ^/'(efc)S anisotropic components with winding num- 
bers ±1 and ±3. Other winding numbers would be of 



higher order in b(kp)T and therefore need not to be con- 
sidered in the diffusive regime. 

Using this expansion of the spin density, the 1st order 
electron-electron collision integral (Eq. (A7|) reduces in 
the diffusive regime to: 



0(1) _ 



= -S X (Xk,lSk,±l +Xk,3Sk,±3), 



(B2) 



where xi an d X3 have been defined in Eq. ( 18 ) 



Plugging the Hartree-Fock precession term (B2) back 



into the equation for the spin density ( 10 ) one finds for 



the isotropic components of the spin density 



v v 

dtSo = -^d x s c --d y s s -s c xb c -s s xh s -s c3 xh c3 ~s s3 xb s3 



with 



(B3) 



S c = Si +S_i, s c3 = s 3 +s_ 3 (B4) 
s s = i(si - s_i), s s3 = i(s 3 - s_ 3 ) 

and the spin-orbit fields 

b c = vpk{— a + P)e y , b s = vpk(a + /3)e x , (B5) 

fc 3 k 3 
b C 3 = -7«F^j-e y , b s3 = jv F — e x . 

In order to obtain a closed equation for So one has to 
determine s c , s s , s c3 and s s3 from the anisotropic compo- 
nents of the Boltzmann equation. In the diffusive regime 
it is sufficient to find the (quasi)-equilibrium solutions for 
the anisotropic coefficients, which are obtained by omit- 
ting the time derivative of the anisotropic components. 
The justification for doing so is that, in order to capture 
the slow precession-diffusion dynamics of the real space 
density, we can interpret the time derivative as a coarse- 
grained one, i.e. d t s — > As /At with At « bp 1 >• r. 
Then the fast relaxation of the anisotropic components 
into the steady state at the beginning of each time in- 
terval At contributes only in higher order in bp t to the 
average over At. Using these approximations one finds 
for the winding number one components 

s c = — Tivd x s + 2Txh c x s + 2tiBi x s c (B6) 
s s = — TivdySo + 2rib s x So + 2tiBi x s s , 



where Bi = \iS is the Hartree-Fock field felt by the 
winding number one spins and T\ is the effective relax- 
ation time for the winding number one part of the spin 
distribution function. Although in our model t\ = r, 
this is not true in general, e.g. consideration of 2nd order 
electron-electron scattering^ will add an additional term 
to T\ . Therefore we use the more general t\ in all our ex- 
pressions in order to make an extension to including e.g. 
2nd order electron-electron interaction straightforward. 



Solving Eqs. (B6) one finds 



s c = -fi {wc^so - 2b c x s + 2ti«(BiX d x s ) (B7) 

-4tiBiX (b c xs ) + 4ti«Bi(Bi • d x s )} 
s s = -fi {vdyS - 2b s x s + 2tiv(BiX d y s ) 



-4riBiX (b s xs ) + 4r 1 2 uB 1 (B 1 • d y s )} 



where we have used Bi||sq and 



n = 



n 



l + (2S 1 r 1 ) 5 



(B8) 



Similarly, the steady state of the winding number three 
anisotropic components of the spin density is governed by 
the following equations 



s c3 
s s3 

T~3 



2b c3 x s + 2B 3 x s c3 



2b s3 x s + 2B 3 x s s3 , 



(B9) 



where B 3 = x 3 S is the Hartree-Fock field felt by the 
winding number three spins and r 3 is the effective relax- 
ation time for the winding number three part of the spin 
distribution function. Although in our model r 3 = r we 
keep the more general r 3 in all our expressions based on 
the same reasoning as detailed for n below Eq. (B6|. 
Solving the equations for s c3 and s s3 one finds 

s C 3 = 2f 3 {b c3 x s + 2r 3 B 3 x (b c3 x s )} (BIO) 
s S 3 = 2f 3 {b s3 x s + 2r 3 B 3 x (b s3 x s )} , 

where we have used that B 3 ||sq and 



T3 = 



T 3 



1 + (2S 3 r 3 ) 



2 ' 



(Bll) 



Plugging the solutions for the anisotropic compo- 
nents (B7) and (BIO I back into the equation for the 
lB3b 



isotropic spin density (|B3j) and using a = /?, i.e., the 
condition for the PSH state, one finds for the spin den- 
sity S, which is given by S = J J^Sq at T = the spin 



diffusion equation Eq. ( 13 ) 
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This reasoning is based on the equations of appendix |b| 
When the expression for the winding number three part of 
the spin distribution (Eq. (B9l) is plugged into the equa- 



tion for the isotropic component (Eq. (B3l) one finds that 
the only term, which generates a finite ^-component is of 
the form b c3 x (s x (b c3 x s )) = -(b e3 x s )(b c3 • s ) 
where b c3 oc e y . 



